Cosmological constraints on DGP braneworld gravity with brane tension 
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We perform a Markov chain Monte Carlo analysis of the self-accelerating and normal branch 
of Dvali-Gabadadze-Porrati braneworld gravity. By adopting a parametrized post-Friedmann de- 
scription of gravity, we utilize all of the cosmic microwave background data, including the largest 
scales, and its correlation with galaxies in addition to the geometrical constraints from supernovae 
distances and the Hubble constant. We find that on both branches brane tension or a cosmological 
constant is required at high significance with no evidence for the unique Dvali-Gabadadze-Porrati 
modifications. The crossover scale must therefore be substantially greater than the Hubble scale 
HoVc > 3 and 3.5 at the 95% C.L. with and without uncertainties from spatial curvature. With 
spatial curvature, the limit from the normal branch is substantially assisted by the galaxy cross 
correlation which highlights its importance in constraining infrared modifications to gravity. 



I. INTRODUCTION 

Cosmological tests of the acceleration of the expan- 
sion offer unique opportunities to test gravity on large 
scales and low curvature. Dvali, Gabadadze, and Por- 
rati (DGP) [l[ proposed that such infrared modifications 
to gravity might arise in a braneworld model where our 
Universe is a 4D brane embedded in a 5D bulk. 

The two branches of cosmological solutions in the DGP 
model have distinct properties. In the so-called self- 
accelerating branch, late-time acceleration of the Uni- 
verse occurs without the need of a cosmological con- 
stant Unfortunately without a cosmological constant, 
the self-accelerating branch predicts cosmological observ- 
ables that are now in substantial conflict with the data 
(e.g. [1, 0, [1]). Moreover, the linearized theory implies 
the presence of ghost degrees of freedom (e.g. The 
former problem can be alleviated with the restoration of 
a cosmological constant or brane tension. A definitive as- 
sessment of the latter problem awaits nonlinear solutions 
[1, On the second or normal branch, self-acceleration 
does not occur but interestingly phantom effective equa- 
tions of state with p/p < —1 can be achieved without 
ghosts with the help of brane tension In both cases, 
brane tension is required but substantial modifications to 
large-scale gravitational dynamics can still persist. 

In this paper, we conduct a Markov chain Monte 
Carlo (MCMC) study of both branches of the DGP 
model using data from cosmic microwave background 
(CMB) anisotropies, supernovae distances, and the Hub- 
ble constant. For observables in the linear regime, we 
adopt the parametrized post-Friedmann (PPF) frame- 
work pH . Il^ and its implementation into a standard 



theoretical predictions. This framework allows us to in- 
clude information from the near horizon scales which are 
crucial for assessing the viability of the self-accelerating 
branch. We also utilize information from the cross corre- 
lation between high-redshift galaxies and the CMB which 
has been proposed as an interesting test of both the self- 
accelerating and normal branches [HI . [l^ . 

In Sec. [ni we review the theory of the normal and 
self- accelerating branches of DGP gravity and their ap- 
proximation through the PPF formalism. We present the 
results of our MCMC study in Sec. IIIII and discuss them 
in Sec. lIVI Finally, the details about the modifications to 
the ISWWLL code [13, [II used for the galaxy-ISW cross- 
correlation observations are specified in the Appendix. 



II. NORMAL AND SELF-ACCELERATING 
BRANCHES 

In the DGP model [l| our Universe is a (3-|-l)-brane 
embedded in a 5D Minkowski space described by the ac- 



tion 
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where 5D quantities are denoted by hats and 4D quan- 
tities are denoted by tildes. Matter fields, including a 
cosmological constant or brane tension and represented 
by Lt, are confined to the brane while only gravity ex- 
tends to the full 5D bulk. We assume that there is no 



Einstein-Boltzmann linear theory solver [l3| for the bulk tension. The constants and /i^ are proportional 



2 



to the inverse Planck masses in the bulk and brane, re- 
spectively. 

Gravity on the brane is consequently modified at large 
scales. In particular, the crossover distance rc ~ 
governs the transition from 5D to 4D scalar-tensor grav- 
ity. On scales smaller than the Vainshtein radius r* = 
(r^Tg)^/'^, nonlinear interactions return gravity to general 
relativity around a point mass with Schwarzschild radius 
Tg. In the following sections we describe the evolution of 
the background and linear density perturbations in the 
two branches of cosmological solutions. 



A. Background expansion 

Variation of the action yields the modified Einstein 
equations on the brane which reduce to the modified 
Friedmann equation in a homogeneous and isotropic met- 



K ^ K 



where H = d/a is the Hubble parameter, K is the spatial 
curvature, a is the scale factor, and pi are the energy 
densities of the various components on the brane. Here 
a — zkl and designates the branch of the cosmological 
solutions. 

For (T = +1, late-time acceleration occurs even with- 
out a cosmological constant A Q and so this branch is 
referred to as self-accelerating DGP (sDGP). In order to 
separate tests of gravity from explanations of accelera- 
tion, we will also study the sDGP branch supplemented 
by a nonvanishing p\ which we denote as sDGP-|-A where 
confusion might arise. For a — —1, DGP modifications 
slow the expansion rate and the branch is referred to as 
normal DGP (nDGP). Here, a cosmological constant is 
required in order to achieve late-time acceleration. 

With the usual definitions ^li ~ p^pi{a = l)/H^ and 
Q,K ~ —K/Hq, the modified Friedmann equation be- 
comes 



H 



where we have assumed that the energy density compo- 
nents include nonrelativistic matter, radiation, and pos- 
sibly a cosmological constant. Here 
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where 
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(5) 



represents the effective contribution of the DGP modifi- 
cation to the energy density assuming the ordinary Fried- 
mann equation. Specifically, 



PDGP = ^ [H^ + §) - Y.P^- 



(6) 



As with any real energy density component, it obeys the 
conservation equation 



PdGP = ^3(1 + t«DGp)PDGP- 

Using Eqs. ([2]) and dZ]), we derive 



1 + WDGP 



(7) 



(8) 



For cases with a cosmological constant it is also useful 
to define the total effective dark energy 



Pe = Pdgp + PA 
and its equation of state 



1 + We = [l + WDGP)- 



Pdgp 



PDGP + PA 



(9) 



(10) 



In nDGP this quantity diverges when the DGP and A 
density terms are equal and opposite at which point the 
value of 1 + We switches sign. In particular, its value 
today is given by 



1 - n 
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(11) 



where we have neglected the small radiation contribu- 
tion. With realistic cosmological parameters wq > —1 
and wq < —1 for sDGP and nDGP, respectively, with 
uiQ = — 1 being the limit of cosmological constant domi- 
nation in either case. 



B. PPF linear theory 

Unlike the background, the evolution of density and 
metric perturbations on the brane require solutions for 
the bulk metric equations. The parametrized post- 
Friedmann framework was introduced in Refs. [13, 
to encapsulate these effects in an effective 3-1-1 descrip- 
tion. Fits to the bulk calculation in sDGP without A or 
curvature from [l^ were given in 11 1 and incorporated 



201 in d]. We extrap- 



into the linear theory code CAMB 
olate these fits to cases with A and curvature here though 
we caution the reader that these have not been tested by 
explicit bulk calculations. For nDGP, we utilize a de- 
scription from [2l| based on bulk calculations from [l^ 
and [l^l with A but no curvature. We again extrapolate 
these results for spatial curvature. The errors induced by 
these extrapolations are controlled given the well-defined 
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limits of A domination and the small dynamical effects 
of curvature in the regime we consider. 

Given the expansion history, the PPF framework is de- 
fined by three functions and one parameter. From these 
quantities, the dynamics are determined by conservation 
of energy and momentum and the Bianchi identities. The 
defining quantities are g(a, k) which quantifies the effec- 
tive anisotropic stress of the modifications and distin- 
guishes the two gravitational potentials, f({a), which de- 
fines the relationship between the matter and the metric 
on superhorizon scales, and fcio-): which defines it in the 
linearized Newtonian regime. The additional parameter 
defines the transition scale between the superhorizon and 
Newtonian behaviors. 

More specifically. 



g{a,k) 



$ + xj/ 



(12) 



where the scalar linear perturbations are represented in 
longitudinal gauge 



-(l + 2^')d^2 + a2^1 + 2$)<ix^ (13) 



where dx^ is the unperturbed spatial line element with 
constant curvature K. In the quasistatic high k limit, 
the DGP model predicts 
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(14) 

where ri_ff(a) = Hn^K/H^a^ ■ On superhorizon scales, 
we take for sDGP 



Hi 



ffSH,sDGPlaj 
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X 1 + 



(15) 



Htc^I - nK{a) - 1.08 



where kn = k/aH, Cg = 0.14 for sDGP and Cg = 0.4 for 
nDGP, respectively. Furthermore, we set Ug ~ 3. 

The function /^(Ina) relates the metric to the density 
at superhorizon scales and is well described by /^(Ina) = 
0.4(7sH(ln a). In the quasistatic regime, the analogous 
relationship between $ — and the density is the Poisson 
equation and that is unmodified from ordinary gravity for 
DGP. Hence fcilna) = 0. 

Finally the parameter cr relates the transition scale in 
the dynamical equations from superhorizon to quasistatic 
behavior. For sDGP we take cr = 1 following [ll[ and 
we employ this value for cases that include A. In nDGP, 
cp ~ 0.15 from [2l| implying a delayed approach to qua- 
sistatic behavior. 



III. CONSTRAINTS ON THE MODELS 

We will use a variety of cosmological data sets 
to constrain the two branches of the DGP mod- 
els. First we use the CMB anisotropy data from 
the five-year Wilkinson Microwave Anisotropy Probe 
(WMAP) [111, the Arcminute Cosmology Bolometer Ar- 
ray Receiver (ACBAR) j25j , the Cosmic Background Im- 
ager (CBI) [21], and the Very SmaU Array (VSA) 
Next we employ data from the Supernovae Legacy Sur- 
vey (SNLS) [2^ and the measurement of the Hubble con- 
stant from the Supernovae and Hq for the Equation of 
State (SHOES) [23| program. Finally we take galaxy- 
ISW (gISW) correlation observations using the likelihood 
code of [U, [l^ ■ We quote results with and without the 
gISW constraint to highlight its impact on the results. 

In Sec. nil Al we discuss the predictions for these ob- 
servables in the two branches of the DGP model. In 
Sees. HITB] and IhTCI we present the resuks of a MCMC- 
likelihood analysis for flat and nonflat universes, respec- 
tively. The MCMC analysis is conducted with the pub- 
licly available CosmoMC [30| package. 



We exclude models — QrHoTc > 1.08 from consider- 
ation as they are not within the domain of applicability 
of the fit nor are they cosmologically viable. For nDGP 
we take [2l[ (c/. (H) 



5SH,nDGp(a) 



1 



2Hrcy/l - nK{a) + 1 



(16) 



The corrections for curvature have not been verified 
by a bulk calculation for the superhorizon cases. For 
the curvatures that we will consider the total impact is 
small as can be verified by omitting the correction. We 
expect therefore that corrections on the correction have 
negligible impact. 

At intermediate scales, g is fitted by the interpolating 
function 



l{a,k) 



1 + (Cgfcff)"' ' 



(17) 



A. Model predictions 

In this section we illustrate model predictions of the 
various cosmological observables we use in the con- 
straints. We chose the parameters of the various models 
that highlight results from the MCMC analysis. 

At high redshifts the DGP modifications become neg- 
ligible on either branch [see Eq. ([3])], and so we choose 
a parametrization that separates high-redshift and low- 
redshift constraints. Specifically we take 6 high-redshift 
parameters: the physical baryon and cold dark matter 
energy density ^l^h'^ and ^Ich'^, the ratio of sound horizon 
to angular diameter distance at recombination multiplied 
by a factor of 100 9, the optical depth to reionization r, 
the scalar tilt n^, and amplitude As at A:* = 0.002 Mpc""'^. 

The low-redshift parameters differ in the various 
classes of models. For flat ACDM and sDGP without A 
there are no additional degrees of freedom. Note that 9 
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carries information on Hq. For flat sDGP+A and nDGP, 
f2A is an extra degree of freedom. For the open versions 
of all models flK is the final degree of freedom. 

For ACDM and sDGP wc illustrate predictions from 
the nonflat maximum likelihood models found in the next 
section (see Tables IVll and |VII| . Since the large-scale be- 
havior of nDGP is new to this work, we highlight the de- 
pendence of observables on Q\ and Qk while keeping the 
high-redshift cosmological parameters fixed (see TablelJ). 
Note in the Tc ^ limit where fi^^ = 0, both nDGP and 
sDGP-l-A become ACDM. We therefore choose to illus- 
trate the maximum likelihood sDGP model with A = 0. 



nDGP 


A 


B 


C 


D 


E 


F 


G 




0.77 


1.00 


1.25 


1.50 


1.50 


1.25 


1.00 


Ok 










-0.025 


-0.015 


-0.010 




0.000 


0.012 


0.049 


0.114 


0.132 


0.057 


0.015 


^^0 


73 


77 


82 


86 


71 


72 


71 



TABLE I: Different choices of nDGP models for illustration. 
Note that nDGP-A is the best-fit (with gISW) flat nDGP 
model, corresponding to ACDM. Other chain parameters are 
fixed to values in Table FVl 



1. Cosmic microwave background 

The CMB probes the geometry of the background ex- 
pansion as well as the formation of large-scale structure. 
The latter manifests itself on the largest scales through 
the integrated Sachs- Wolfe (ISW) effect from the evolu- 
tion of the gravitational potential. To predict these ef- 
fects we implement the PPF modifications from Sec. Ill Bl 
The incorporation of the PPF formalism into a standard 
Einstein-Boltzmann linear theory solver yields an effi- 
cient way to obtain predictions of the DGP model for the 
CMB. We utilize the PPF modifications to CAMB ^ 
implemented in Refs. [l^l, which we can apply di- 
rectly for sDGP and figure as a starting point for the 
implementation of nDGP and sDGP-|-A. In Fig. [1] we 
plot the CMB temperature anisotropy power spectrum 
with respect to angular multipole £ for the best-fit mod- 
els of ACDM and sDGP, as well as the nDGP parameter 
choices given in Table [D 

Relative to ACDM, the growth of structure is sup- 
pressed in the sDGP model, yielding an ISW enhance- 
ment at the lowest multipoles. This enhancement is suf- 
ficiently large to bring the sDGP model without A into 
serious conflict with the joint data Q. The opposite ef- 
fects occurs in the nDGP model and lead to predictions 
that are compatible with CMB data. Here raising J^a at 
fixed Hk enhances the low multipoles through the ISW 
effect. However, compensating effects from curvature can 
lead to parameter degeneracies. 

At high redshifts the contribution of ilr^ to the Hubble 
parameter, Eq. ([3]), becomes negligible in either branch. 
The CMB acoustic peaks can therefore be utilized as 



usual to infer constraints on the high-redshift parameters, 
in particular, the physical energy densities of baryonic 
matter and dark matter as well as the angular diameter 
distance to recombination. 



2. Distances to the supemovae and Hq 

The comparison of the magnitudes of high-redshift to 
low-redshift supernovae yields a relative distance mea- 
sure. Theoretical predictions for the distance modulus 
are related to the luminosity distance, dL{z) = (1 -I- 
z)r(z), where r{z) is the comoving angular diameter dis- 
tance defined by 

f sin [HoV^xiz)] /Ho^WhTl, < 0, 
= \ X{z), = 0, 

[ sinh [HoV^xiz)] /Ho^W^l Hk > 0, 

(18) 

where the comoving radial distance x is 

r dz' 

The supernovae magnitudes, once standardized, are re- 
lated to the distance by 

TO = /i + i\/ = 51ogiorfL+^/ + 25, (20) 

where is in units of Mpc. The unknown absolute 
magnitude M of the supernovae is a nuisance parameter 
in the fit and is degenerate with Hq. Hence supernovae 
measure relative distances within the set. 

In Fig. [21 we plot the predictions for the distance 
modulus for the SNLS data in sDGP gravity, nDGP-B, 
nDGP-F, and in the ACDM model. 

The acoustic peaks in the CMB and the measurement 
of the local Hubble constant additionally provide ab- 
solute distance probes which complement the relative 
distance measure of the supernovae. For the Hubble 
constant, we take the SHOES measurement of Hq = 
74.2 ± 3.6 km s~^ Mpc~^ which employs Cepheid mea- 
surements to link the low-redshift supernovae to the dis- 
tance scale established by the maser galaxy NGC 4258. 



3. Galaxy-ISW cross correlations 

The correlation between galaxy number densities and 
the CMB anisotropics can be used to isolate the ISW ef- 
fect in the CMB. The enhanced ISW effect exhibited in 
the sDGP model without A leaves a strong imprint on 
the large scales of the CMB temperature anisotropy. As 
pointed out by Song et al. , an interesting consequence 
of this is a considerable correlation of high-redshift galax- 
ies with the CMB. 

For nDGP gravity, whereas the ISW effect does not 
exhibit a substantial impact on the largest scales in the 



5 





z 



FIG. 2: Best-fit distance modulus for sDGP, as well as tiie 
overlapping predictions for nDGP-B and nDGP-F with re- 
spect to ACDM. 



CMB, useful signatures remain in the correlations with 
galaxies that can break parameter degeneracies (l6j . 

We evaluate the gISW cross correlations in the Lim- 
ber and quasistatic approximat ion , as it is done in the 
ISWWLL code [13, [ll used for the data analysis. There- 
fore, we write 

= '";f;g;r /-i--/.(--)g(')p(') 

Here, D{z) is the linear growth rate in the quasistatic 
regime defined by A„(fc,z) = Amik,0)D{z)/ D{0), 
where Am(fc, z) is the matter density perturbation. P{k) 
is the matter power spectrum today. 

The approximations in Eq. (j2ip become accurate at 
the percent level for £ > 10. This condition is satisfied by 
about 90% of the total 42 data points that are used in the 



ISWWLL code. We discuss details about the approxima- 
tions and the data in the Appendix. The data are divided 
into nine galaxy sample bins j, i.e., 2MASS0-3, LRGO-1, 
QSOO-1, and NVSS. The function relates the mat- 

ter density to the observed projected galaxy overdensity 
with fj{z) = bj{z)Ilj{z) in the absence of magnification 
bias. Ilj{z) is the redshift distribution of the galaxies 
and the bias factor bj{z) is assumed independent of scale, 
but dependent on redshift. The code determines fj{z), 
among other things, from fitting autopower spectra and 
cross-power spectra between the samples. 

We modify the above calculations in the ISWWLL 
code with the appropriate DGP quantities such that the 
correct predictions for the crosscorrelations are obtained. 
We refer to the Appendix for details. The predictions for 
the best-fit values, combining all data, of ACDM and 
sDGP for the different samples are shown in Fig. [31 We 
also plot the curves for nDGP-B and nDGP-F to illus- 
trate the breaking of the degeneracy between Q\ and 
^Ik- Notice that the model with larger curvature has 
reduced correlation especially at high redshift. We shall 
see that models with significantly larger curvature can be 
excluded by the gISW data. 



B. Flat universe constraints 

We begin by studying a universe without spatial cur- 
vature, where the basic cosmological parameter set is 
P = {Qbh'^, 0, r. Us, ln[10^°As]}. We implement the 

following flat priors on them: rif,/i^ € (0.01,0.1), fic/i^ G 
(0.05,0.99), e e (0.5,10), r e (0.01,0.8), n, € (0.5,1.5), 
and ln[10i°As] G (2.7,4). For nDGP and sDGP with 
nonvanishing A, we use f^A G (0.0,2.5). 

We begin with the analysis of flat ACDM without DGP 
modifications in Table [Til We show constraints with and 
without the gISW data and the maximum likelihood pa- 
rameters and value. Horizontal lines divide the chain 
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FIG. 3: Best-fit ACDM and sDGP galaxy-ISW cross correlations for the different galaxy samples, roughly ordered in increasing 
effective, bias-weighted, redshift. Note the distinct predictions for the previously degenerate nDGP-B and nDGP-F models. 



Parameters 


ACDM 


ACDM (with gISW) 


fOOfift/i^ 


2.248 ± 0.055 


2.240 


2.251 ±0.055 


2.258 




0.1080 ± 0.0043 


0.1072 


0.1075 ±0.0042 


0.1071 


6 


1.0410 ± 0.0027 


1.0404 


1.0411 ±0.0027 


1.0417 


T 


0.086 ±0.017 


0.086 


0.087 ±0.017 


0.089 


Us 


0.963 ±0.013 


0.961 


0.963 ±0.013 


0.965 




3.176 ±0.041 


3.177 


3.174 ±0.041 


3.173 




0.751 ±0.019 


0.754 


0.754 ±0.019 


0.758 




0.249 ±0.019 


0.246 


0.246 ±0.019 


0.242 


Ho 


72.6 ± 1.8 


72.6 


72.8 ± 1.8 


73.2 


-21nL 


2834.29 


2867.99 



TABLE II: Means, standard deviations (left subdivision of 
columns), and best-fit values (right subdivision of columns) 
with likelihood for the fiat ACDM model using data from 
WMAP, ACBAR, CBI, VSA, SNLS, and SHOES without (left 
column) and with the gISW data (right column). 



parameters from the derived parameters and the best-fit 
(maximum) likehhood. In the case of ACDM, the inclu- 
sion of the giSW data does not yield noticeable improve- 
ment on the parameter constraints [l7| . This analysis 
sets the baseline by which adding the DGP degrees of 



Parameters 


sDGP 


sDGP (with g 


ISW) 




2.390 ±0.066 


2.393 


2.390 ± 0.065 


2.376 




0.0884 ± 0.0042 


0.0873 


0.0889 ± 0.0041 


0.0899 


6 


1.0448 ±0.0028 


1.0447 


1.0449 ± 0.0028 


1.0452 


T 


0.105 ±0.021 


0.110 


0.105 ±0.021 


0.103 




1.011 ±0.015 


1.013 


1.011 ±0.015 


1.007 




3.001 ±0.045 


2.998 


3.003 ± 0.044 


3.015 




0.1410 ±0.0075 


0.1430 


0.1403 ± 0.0075 


0.1384 




0.249 ±0.020 


0.244 


0.251 ±0.020 


0.256 


Ho 


67.2 ± 1.7 


67.6 


67.1 ± 1.7 


66.7 


-2AlnL 


32.70 


33.06 



TABLE III: Same as Table HIl but for the fiat sDGP model. 
— 2AlnL is quoted with respect to the maximum likelihood 
fiat ACDM model. 



freedom should be measured. 

In the flat sDGP model without A, there is no choice 
of parameters that can satisfy the joint requirements of 
geometrical measurements from the CMB, supernovae, 
and Hq and the dynamical requirements from the ISW 
effect. For sDGP, we find — 2Alni = 32.7 with respect 
to ACDM and ~2AlnL = 33.1 {5.8a) when including the 
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Parameters 


sDGP+A 


sDGP±A (with gISW) 




2.265 ± 0.058 


2.245 


2.265 ±0.058 


2.257 




0.1050 ± 0.0046 


0.1071 


0.1048 ±0.0046 


0.1070 


e 


1.0415 ± 0.0028 


1.0405 


1.0415 ± 0.0027 


1.0415 


T 


0.089 ±0.017 


0.080 


0.089 ±0.017 


0.083 


TLs 


0.969 ± 0.014 


0.961 


0.969 ±0.014 


0.968 


\n[10^^As] 


3.153 ± 0.044 


3.165 


3.152 ±0.044 


3.154 


O A 


0.590 - 0.752 


0.733 


0.588 - 0.751 


0.719 




< 0.0178 


0.0001 


< 0.0186 


0.0003 




0.248 ±0.019 


0.248 


0.247 ±0.018 


0.247 


Ho 


71.9 ± 1.9 


72.3 


71.9 ± 1.9 


72.5 


-2AlnL 


0.20 


0.13 



TABLE IV: Same as Table HIl but for the flat sDGP±A 
model, quoting one-sided ID marginahzed upper 95% CL for 
fir, and 68% MCI for Qa- 



Parameters 


nDGP 


nDGP (with gISW) 


lOOfift/i^ 


2.237 ± 0.054 


2.245 


2.238 ± 0.056 


2.254 




0.1109 ± 0.0049 


0.1095 


0.1100 ±0.0046 


0.1076 


6 


1.0406 ± 0.0027 


1.0410 


1.0407 ±0.0027 


1.0409 


T 


0.084 ± 0.016 


0.084 


0.085 ±0.017 


0.092 




0.958 ± 0.012 


0.961 


0.959 ±0.013 


0.961 




3.196 ± 0.043 


3.182 


3.191 ±0.043 


3.190 




0.754 - 0.934 


0.765 


0.753 - 0.924 


0.772 




< 0.0228 


0.0001 


< 0.0203 


0.0001 




0.247 ± 0.019 


0.253 


0.243 ±0.018 


0.244 


Ho 


73.6 ± 2.0 


72.2 


73.9 ± 2.0 


73.0 


-2 A hiL 


0.05 


0.23 



TABLE V: Same as Table HVl but for the flat nDGP model. 



gISW likelihood. In this case, the ISW effect is so large 
at low multipolcs that the CMB alone rules out such con- 
tributions Q and the gISW constraint adds only an in- 
significant amount of extra information (see Tables |TTT|. 
The strengthening of the constraint when compared to 
Ref. 0] comes from the improved Hubble constant mea- 
surements. 

In the sDGP+A and nDGP models, the cosmological 
constant becomes a free parameter and we have to add it 
to the parameter set, hence, P ^ P U {Oa}- is a de- 
rived parameter and in particular we get fir, — s- in the 
limit f^A ^ (1 ~ ^m)- In this limit, the phenomenology 
of ACDM is recovered for all observables. Preference for 
a finite fir, indicates evidence for the DGP modification 
in these cases. 

In both the nDGP and sDGP+A cases the maximum 
likelihood models differ insignificantly from ACDM (see 
Tables IIVI and |V| and there is no preference for finite 
fir,- Conversely, both branches require a finite Q\ at 
high significance. 

Since ACDM is the fir, limit of both branches 
with A, the slightly poorer fit for nDGP and sDGP+A 
should be attributed to sampling error in the MCMC. 
The one-sided ID marginalized upper 95% confidence 



limits for fir, are fir, < 0.0178(0.0186) for sDGP±A 
and fir, < 0.0228(0.0203) for nDGP where the values in 
parentheses include the gISW constraint. These values 
indicate that the crossover scale is at least substantially 
greater than the Hubble scale Horc ^ 3.5. 

In this ACDM hmit, the modifications to the gISW 
predictions do not affect the constraints. The slight 
weakening of the constraints with the inclusion of gISW 
in sDGP±A does not indicate a statistically significant 
tension but does suggest that future improvement in con- 
straints can tighten the bounds on HqTc- In particular, 
sDGP modifications tend to enhance correlations at high 
rcdshift relative to low redshift. The current data have 
a marginal preference for increased correlation with red- 
shift relative to ACDM (see Fig. [3]). 

Note that due to the distinctive skewness of the poste- 
rior distribution, we give the ID marginalized 68% mini- 
mum credible intervals (MCI) (see Ref. [3l| ) for the brane 
tension fl{^ as opposed to the standard deviations given 
for the other parameters. 

Finally, in the context of these fiat models the possi- 
bility of phantom equations of state currently is highly 
constrained. For nDGP I + wq > -0.039 at the 95% CL. 



C. Nonflat universe constraints 

In a universe with spatial curvature, we include fix as 
a parameter in the chain for each of the model classes. 
We use the prior flK & (—0.1,0.1), which we weaken to 
fix G (—1, 1) in nDGP since we expect degeneracies be- 
tween fix and flA. We also implement latter prior for 
sDGP±A. For ACDM, Ho et al. ^ have found an im- 
provement of the constraints on fix by a factor of 3.2, 
with respect to WMAP3 data alone, due to the inclusion 
of the gISW and weak Icnsing data. However we find 
that the inclusion of the other data, specifically the su- 
pernova and Hq data, make curvature constraints only 
marginally improved by the gISW inclusion. We again 
use these ACDM results shown in Table I VII as a base- 
line for comparison with sDGP, sDGP+A, and nDGP in 
Tabslcs Eni IVllIl and |TXl 

For sDGP without A, adding curvature alleviates the 
tension between CMB and supernova distance measures. 
However, it cannot reduce the ISW contributions [T3 | 
and so we obtain — 2AlnL = 23.3(23.8), with respect 
to ACDM where values in parentheses include the gISW 
constraint. Utilizing all of the data, the significance of 
the exclusion of sDGP without A is ~ 5a. 

Similar to the flat case, we find no preference for a 
finite fir, in nDGP and sDGP+A and consequently no 
indications of DGP modifications to gravity (see Fig.d]). 
With sDGP±A, we are again driven to the limiting case 
of ACDM with the slightly poorer best fit reflecting a 
sampling error in the chain. Allowance for curvature on 
the other hand weakens the upper limit on the DGP mod- 
ifications: fir, < 0.0248(0.0244) and HoTc > 3.18(3.20) 
at 95% CL. 
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For nDGP; the addition of curvature introduces a de- 
generacy with the cosmological constant. As was pointed 
out by Giannantonio et ah this degeneracy can be 
broken by the use of galaxy-ISW cross correlations since 
high curvature sohitions underpredict the correlation es- 
pecially at high rcdshift. Figure [5] illustrates this de- 
generacy and the effect of gISW measures. The result 
of marginalizing curvature in nDGP is again a weaken- 
ing of the DGP constraints Qr, < 0.0501(0.0300) and 
Hor^ > 2.23(2.89) at 95% C.L. ' 

In summary with the gISW constraint, the limit on 
either branch implies HqTc ^ 3 and only a small weak- 
ening from the flat case of 3.5. Furthermore due to the 
curvature degeneracy in nDGP, restrictions on phantom- 
like equations of state are also somewhat weakened to 
wo + 1 < -0.049. 



Parameters 


ACDM 


ACDM (with gISW) 




2.250 ±0.056 


2.246 


2.249 ± 0.055 


2.238 




0.1084 ±0.0052 


0.1095 


0.1084 ± 0.0051 


0.1085 


e 


1.0412 ± 0.0027 


1.0412 


1.0411 ± 0.0027 


1.0419 


T 


0.086 ±0.017 


0.090 


0.087 ± 0.017 


0.083 


n. 


0.963 ±0.013 


0.960 


0.963 ± 0.013 


0.962 


ln[10"^.] 


3.176 ±0.044 


3.196 


3.179 ± 0.043 


3.174 




-0.0001 ±0.0063 


0.0020 


0.0007 ± 0.0062 


0.0021 




0.751 ±0.020 


0.751 


0.753 ± 0.019 


0.758 




0.249 ± 0.022 


0.246 


0.246 ± 0.022 


0.240 


Ho 


72.6 ±3.0 


73.2 


73.0 ±3.0 


73.8 


-2\nL 


2834.01 


2867.74 



TABLE VI: ACDM as in Table HI except allowing spatial 
curvature. 



Parameters 


sDGP 


sDGP (with g 


ISW) 




2.377 ± 0.061 


2.365 


2.376 ±0.062 


2.352 




0.0951 ± 0.0041 


0.0970 


0.0952 ±0.0039 


0.0979 


6 


1.0441 ± 0.0028 


1.0451 


1.0441 ±0.0028 


1.0439 


T 


0.091 ± 0.020 


0.084 


0.092 ±0.019 


0.084 




1.004 ±0.014 


1.002 


1.004 ±0.014 


0.997 


ln[10^°A.] 


3.018 ± 0.043 


3.019 


3.021 ±0.043 


3.037 




0.0186 ± 0.0055 


0.0212 


0.0182 ±0.0055 


0.0220 




0.1486 ± 0.0068 


0.1486 


0.1479 ±0.0067 


0.1467 




0.218 ±0.019 


0.216 


0.220 ±0.019 


0.220 


Ho 


74.0 ± 3.0 


74.7 


73.7 ±2.9 


74.2 


-2AlnL 


23.32 


23.79 



TABLE VIL sDGP without A as in Table|ITIl except allowing 
spatial curvature. — 2AlnL is quoted with respect to the 
maximum likelihood ACDM model with curvature here and 
in the following tables. 



IV. DISCUSSION 

We have performed the first Markov chain Monte 
Carlo analysis of the nDGP and sDGP branches of DGP 



Parameters 


sDGP±A 


sDGP±A (with 


gISW) 


lOOSlfc/i^ 


2.266 ±0.058 


2.252 


2.266 ± 0.0059 


2.251 




0.1065 ±0.0051 


0.1066 


0.1064 ± 0.0051 


0.1095 


e 


1.0416 ±0.0028 


1.0406 


1.0415 ± 0.0028 


1.0414 


T 


0.087 ±0.017 


0.077 


0.088 ± 0.017 


0.090 


Tis 


0.968 ±0.014 


0.962 


0.968 ± 0.014 


0.960 


ln[10^"As] 


3.157 ±0.046 


3.153 


3.158 ± 0.046 


3.197 




0.0032 ± 0.0068 


0.0022 


0.0036 ± 0.0065 


0.0018 


O A 


0.557 - 0.745 


0.711 


0.561 - 0.746 


0.737 




< 0.0248 


0.0006 


< 0.0244 


0.0000 




0.245 ±0.022 


0.240 


0.243 ±0.021 


0.248 


Ho 


72.8 ± 3.0 


73.3 


73.1 ±2.9 


72.9 


-2AlnL 


0.07 


0.04 



TABLE VIIL sDGP with A as in Table HVl but allowing spa- 
tial curvature. 



Parameters 


iiDGP 


nDGP (with gISW) 




2.239 ±0.056 


2.245 


2.242 ±0.056 


2.239 




0.1099 ±0.0054 


0.1076 


0.1094 ±0.0054 


0.1099 


e 


1.0409 ± 0.0027 


1.0412 


1.0409 ±0.0027 


1.0409 


T 


0.084 ±0.017 


0.084 


0.085 ±0.017 


0.091 


Us 


0.959 ±0.013 


0.960 


0.960 ±0.013 


0.960 




3.189 ±0.045 


3.176 


3.188 ±0.045 


3.205 




-0.0055 ± 0.0080 


-0.0056 


-0.0029 ± 0.0069 


0.0021 


Ha 


0.749 - 1.009 


0.801 


0.749 - 0.953 


0.764 




< 0.0501 


0.0008 


< 0.0300 


0.0000 




0.255 ±0.023 


0.261 


0.248 ± 0.022 


0.247 


Ho 


72.1 ±3.0 


70.6 


73.0 ± 3.0 


73.1 


-2AlnL 


0.09 


0.41 



TABLE IX: nDGP as in Table |Vl but allowing spatial curva- 
ture. 



braneworld gravity to utilize all of the CMB data, in- 
cluding the lowest multipoles, and its correlation with 
galaxies (gISW). We also include supcrnovae and Hub- 
ble constant data in the constraint. 

We find no preference for DGP modifications to gravity 
on either branch. Indeed, on the self-accelerating branch 
without A, the model is excluded at the 4.9(t and S.Scr 
levels with and without curvature respectively . While 
the gISW data do not substantially improve this con- 
straint, they do additionally disfavor sDGP. 

With the inclusion of A on either branch, the DGP 
model cannot be entirely excluded but its modifications 
are strongly limited. We find that the crossover scale, 
which measures the strength of the modifications, must 
be substantially above the Hubble scale H^rc > 3 with 
curvature and 3.5 without curvature. The robustness of 
this constraint is substantially assisted by the gISW data. 
In nDGP, it breaks the geometric degeneracy between A 
and spatial curvature. In sDGP, the relatively large cor- 
relation at high redshift offers opportunities in the future 
for improving the limits on HqTc- These abilities high- 
light the importance of obtaining improved gISW data 
for constraining infrared modifications to gravity. 
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FIG. 4: Marginalized likelihood for Slr^ in the nonflat nDGP and sDGP models. The vertical lines indicate 68%, 95%, and 
99% C.L. 




FIG. 5; Contours of 2D marginalized 68%, 95%, and 99% confidence boundaries using WMAP5, ACBAR, CBI, VSA, SNLS, 
and SHOES (left panel), including gISW (right panel) for nDGP. 
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We use the publicly available ISWWLL code [l7|, llSf 
for our analysis. Note that we have turned off weak lens- 
ing likelihood contributions in the code, focusing only 
on the gISW constraints. The 42 data points of giSW 
cross correlations that are used in the likelihood analy- 
sis are collected from the Two Micron All Sky Survey 
(2MASS) extended source catalog (XSC) 0, [H, the 
luminous red galaxies (LRG) and photometric quasars 
(QSO) of the Sloan Digital Sky Survey (SDSS) [H, 
and the National Radio Astronomy Observator y (N RAO) 
Very Large Array (VLA) Sky Survey (NVSS) [3|. They 
are divided into nine galaxy sample bins j (2MASS0- 
3, LRGO-1, QSOO-1, NVSS) based on flux (2MASS) or 
redshift (LRG, QSO). These data points are a selection 
of multipole bins from all samples, where the selection 
is based on the avoidance of nonlincarities and system- 
atic effects from dust extinction, galaxy foregrounds, the 
thermal Sunyaev-Zel'dovich effect, and point source con- 



10 



tamination to affect the gISW cross correlations [l7l |. 

In the remainder of this Appendix, we discuss the de- 
tails of the modifications implemented in the ISWWLL 
code. First, we describe the calculation of the quasistatic 
linear growth rate D{z) in the giSW cross correlation, 
Eq. ([21]) . and analyze the validity of the Limber and the 
quasistatic approximation. We then discuss the function 
fj{z) that carries information about the redshift distri- 
bution and bias. 



1. gISW cross correlations 

It has been argued that for nDGP and sDGP the 
gISW cross correlations are well described within the 
quasistatic regime [13, [H, [l^, HI]. Here, this can eas- 
ily be seen from the substitution k {£ + l/2)/x(z) 
considering the relevant redshifts. In this limit, we solve 
the ordinary differential equation [s^] 

A:'„+ (^2 + A'^- ^ (1 - .gQs) f|§?A,„ = (Al) 

for the linear matter density perturbation . Note that 
for nDGP, in the limit Vc oo, we have gqg and 
H{z) approaches the expansion history of ACDM. There- 
fore, in this limit, Eq. (jAip recovers the quasistatic or- 
dinary differential equation for the matter overdensity 
in ACDM. We solve Eq. (|Aip with initial conditions at 
ai <C 1, in a regime where general relativity is expected 
to hold, i.e., A^(ai) = Am{ai) with a normalization set 
by the initial power spectrum. 

The accuracy of the Limber approximation in the case 
of ACDM is at the order of 10% at € = 2 and drops 
approximately as i"^ at higher £ (see e.g. [H, [H, |4C|). 
The error depends further on the width of the redshift 
distribution, which changes only little with DGP effects. 
The relative deviation from the exact result at ^ = 6 
docs not exceed ^ 3% for the samples and typical models 
considered in Fig. [3] Given the large errors of the cur- 
rently available data points at low i, we conclude that 
the Limber approximation is applicable and furthermore 
very useful since it is numerically faster than an exact 
integration. 



the Markov chain, fj{z) is recomputed when changing the 
cosmological parameter values. The methods by which 
this function is determined differ for each sample, but 
they are all based on galaxy clustering data. 

The 2MASS galaxies are matched with SDSS galax- 
ies in order to identify their redshifts. To obtain the 
nonlinear power spectrum, the Q model for nonlincari- 
ties [4l| is applied. Then, the code computes the galaxy 
power spectrum and fits it to measurements, thereby de- 
termining the bias b{z) and Q. Since the required ac- 
curacy for the estimation of bias is only at the order 
of a few tens of percent [l3|, this processing is also ap- 
plicable to DGP. The Q model is also adopted for LRG 
galaxies, where the redshift probability distribution is in- 
ferred with methods described in Rcf. [4^. For QSO, 
first, a preliminary estimate for the redshift distribution 
is deduced by locating a region of sky with high spec- 
troscopic completeness, but simultaneously maintaining 
a large area. Taking into account magnification bias and 
fitting bj{z)Ilj(z) using the quasar power spectrum and 
quasar-LRG cross power yields the desired shape of fj (z) . 
Finally, the effective redshift distribution of NVSS is ob- 
tained from cross-correlating with the other samples and 
fj{z) is fitted with a P distribution. 

The part of the ISWWLL code that is devoted to this 
processing is configured for a parametrization of the ex- 
pansion history by We = wq + {I — a)wa- This approach 
gives a good approximation to sDGP in the domain of 
interest, but it fails for nDGP and sDGP+A due to the 
appearance of a divergence in We{a). Therefore, instead 
of taking wq and Wa to describe the expansion history, 
we utilize fi^c a-nd f2A, where only the latter really is a 
necessary, free parameter. 

In case of the SDSS quasars, the derivation of fj{z) in- 
volves the linear growth factor, which we need to replace 
by its DGP counterpart. This implies solving Eq. (jAl[) . 
A further contribution for the QSO samples is due to 
magnification bias. In the quasistatic regime of DGP the 
relationship between the metric combination sensitive to 
gravitational redshifts and lensing ($ — and the den- 
sity perturbations is unmodified so the expression of the 
lensing window function for magnification effects given 
in Ref. [l3| is unchanged. 



2. Redshift distribution and bias 



A further modification to the code that we need to 
conduct is in the determination of the function fj{z). In 
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